Transport and Percolation Theory in Weighted Networks 
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We study the distribution P(a) of the equivalent conductance a for Erdos-Renyi (ER) and scale- 
free (SF) weighted resistor networks with N nodes. Each link has conductance g = e~ ax , where 
x is a random number taken from a uniform distribution between and 1 and the parameter a 
represents the strength of the disorder. We provide an iterative fast algorithm to obtain P(a) and 
compare it with the traditional algorithm of solving Kirchhoff equations. We find, both analytically 
and numerically, that P(a) for ER networks exhibits two regimes, (i) A low conductance regime 
for a < e~ ap " where p c — is the critical percolation threshold of the network and (k) is 

average degree of the network. In this regime P(cr) is independent of N and follows the power law 
P(cr) ~ o~ a , where a = 1 — (k)/a. (ii) A high conductance regime for a > e~ aVc in which we find 
that P(o~) has strong N dependence and scales as P(er) ~ f(a, apc/N 1 ^ 3 ). For SF networks with 
degree distribution P(k) ~ k~ x , k m i n < k < k max , we find numerically also two regimes, similar to 
those found for ER networks. 



Recently much attention has been focused on complex 
networks which characterize many biological, social, and 
communication systems [E S S 3- The networks are 
represented by nodes associated with individuals, orga- 
nizations, or computers and by links representing their 
interactions. In many real networks, each link has an as- 
sociated weight, the larger the weight, the harder it is to 
transverse the link. These networks are called "weighted" 
networks [5|, [y] . 

Transport is one of the main functions of networks. 
While the transport on unweighted networks has been 
studied p} , the effect of disorder on transport in networks 
is still an open question. Here we study the distribution 
P(cr) of the equivalent electrical conductance a between 
two randomly selected nodes A and B on Erdos-Renyi 
(ER) @,[|[ and scale-free (SF) [l[ weighted networks. We 
first provide an iterative fast algorithm to obtain P(a) 
for disordered resistor networks, and then we develop a 
theory to explain the behavior of P(a). The theory is 
based on the percolation theory [10( for a weighted ran- 
dom network. We model a weighted network by assigning 
the conductance of a link connecting node i and node j 
as in Ref. [Tfll 
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where the parameter a controls the broadness 
("strength") of the disorder, and random 
number taken from a uniform distribution in the range 
[0,1]. We use this kind of disorder since a recent study 
of magnetorresistance in real granular materials systems 
[IH shows that the conductance is given by Eq. |T]). 
Moreover, a recent study [HJ shows that many types of 
disorder distributions lead to the same universal behav- 
ior. The range of a 3> 1 is called the strong disorder 
(SD) limit [lj, The special case of unweighted 

networks, i.e., a = or gij = 1 for all links have been 
studied earlier [7[. 



To construct ER networks of size N, we randomly con- 
nect nodes with (k)N/2 links, where (k) is the average 
degree of the network. To construct SF networks, in 
which the degree distribution follows a power law, we 
employ the Molloy-Reed algorithm The traditional 
algorithm to calculate the probability density function 
(pdf) P{a) is to compute a between two nodes A and B 
by solving the Kirchhoff equations with fixed potential 
Va = 1 and Vb = and compute P(a)da, which gives 
the probability that two nodes in the network have con- 
ductance between a and a + da. However, this method 
is time consuming and limited to relatively small net- 
works. Here we also use an iteration algorithm proposed 
by Grimmett and Kesten [l6| to calculate P(a) and show 
that it gives the same results as the traditional Kirchhoff 
method. 
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FIG. 1: Schematic Iteration model. In this example R\ is 
infinite, so it is not taken into account in the sum in Ri of 
Eq. ©. 



In the limit N — > oo we ignore the loops between 2 ran- 
domly chosen nodes because the probability to have loops 
is very small. Hence the resistivity Ri of a randomly se- 
lected branch i connecting a node with infinitely distant 

nodes satisfies Ri = ri + RJ 1 )' wnere r % — e<1Xi 

is the random resistance of the link outgoing from this 
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node and fc is a random number taken from the distri- 
bution pfc = Pk ■ k/(k), which is the probability that a 
randomly selected link ends in a node of degree k, where 
Pk is the original degree distribution. In Fig. [TJ we show 
the schematic iteration method. The randomly selected 
nodes A and B are connected to the infinitely distant 
nodes C. When we calculate Rac, the resistance between 
A and C, we perform the iterative steps as follows: 

First we calculate the distribution of resistivities of the 
branches connecting node A with C. We start with M 
branches having resistivities i?- ^ =0 (i = 1, 2, AT), 
where Af is an arbitrary large number. Thus, initially 
the histogram of these resistivities Pq(R) = 8(R). At the 
iterative step n+1, we compute a new histogram P n+ i(R) 
knowing the histogram P n (R). In order to do this we 

generate a new set of resistivities by connecting in 

parallel fc — 1 outgoing branches coming from a randomly 
selected node of degree fc obtained from the distribution 
Pk = Pk ■ k/(k). Then we compute the resistivity of a 
branch going through this node via an incoming link with 
a random resistivity taken from the link resistivity 
distribution, 
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(n+l) 
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In Eq. ([3]), if at least one of the terms K[ — 0, we take 
R {n+i) = r {n)_ Thug after the firgt iterative gtep Pl (ty 

coincides with the distribution of link resistivities. 

According to the theorem proved in [l6j . as n — > oo, 
P n (R) converges to a distribution of the resistivities of a 
branch connecting a node to the infinitely distant nodes. 
The resistivity between a randomly selected node of de- 
gree k and the infinitely distant nodes is defined by 



R 



(0 



(3) 



where k is selected from the original degree distribution 
Pk and Rj is selected from P n ^oo(R)- 

Finally, to compute the resistivity Rij between two 
randomly selected nodes i and j (for example A and B 
in Fig. [T]), we compute R^ = R^ + R(j), where R^ and 
are randomly selected resistivities between a node 
and the infinitely distant nodes. If N is a sufficiently 
large number, we find the conductance distribution P(o~) 
between any two randomly selected nodes. 

In Figs.^a) and[^b) we show the results of P{u) using 
the traditional method of solving KirchhofFs equations 
for different values of A and the iterative method with 
A — > oo for both ER and SF networks. We see that the 
main part of the distribution (low conductances) does 
not depend on A, and only the high conductance tail 
depends on A. 

The behavior of the two regimes, low conductance and 
high conductance, can be understood qualitatively as 
follows: For strong disorder a ^ 1 all the current be- 
tween two nodes follows the optimal path between them. 
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FIG. 2: Plots of P(u) for several values of N. The symbols are 
for the KirchhofF method and the solid line is for the iterative 
method with N — > oo. (a) ER networks with fixed (fc) = 3 
and a — 15. (b) SF networks with fixed A = 3.5, fc m i n = 2, 
(fc) ~ 3.33 and a — 20. The dashed line slopes are from the 
prediction of Eg. (111) or ([T3J. 



The problem of the optimal path in a random graph in 
the strong disorder limit can be mapped onto a perco- 
lation problem on a Cayley tree with a degree distribu- 
tion identical to the random graph and with a fraction p 
of its edges conducting [l7|. However, the conductance 
on this path is determined by the bond with the lowest 
conductance e~ aXmax , where x max is the maximum ran- 
dom number along the path. In the majority of cases 
Xmax > Pa where p c is the critical percolation threshold 
of the network, and only when the two nodes both belong 
to the incipient infinite percolation cluster (IIPC) [l(|, 
Xmax < Pc- Since the size of the IIPC scales as A 2 / 3 , the 
probability of randomly selecting a node inside the IIPC 
is proportional to N 2 / 3 /N = A" 1 / 3 Then the 

probability of randomly selecting a pair inside the IIPC 
is proportional to (A -1 / 3 ) 2 = A~ 2 / 3 . These nodes con- 
tribute to the high conductance range a > e~ aPa of P{cr). 
The low conductance regime is determined by the distri- 



bution Of Xr> 



that follows the behavior of the order 



parameter P^ (p) (for p > p ^) in the percolation problem 
which is independent of A [Tj]. (This will be explained 
later in the theoretical approach for the low conductance 
regime.) 

We call the low conductance regime a non-percolation 
regime and the high conductance regime a percolation 
regime. In contrast, the property of existing two regimes 
does not show up in the optimal path length [l8l . Il9j and 
only the scaling regime with A appears. This is since 
the path length for almost all pairs is dominated by the 

iipc 0. 

In Figs. [31(a) and[3Jb) we plot for a given A only the 
non-percolation part of P(<r) as a function of a for fixed 
values of (k)/a and different (fc) and a values for ER 
networks. We see that it obeys a power law with the 
slope (k)/a — 1 for a < e~ ap ". Note that for ER networks 
Pc = l/(k) [1, Q. In Fig. |3](c), we plot the conductance 
distribution for SF networks for fixed values of (k)/a. We 
can see the non-percolation part seems to obey the same 
power law as ER networks. 

Next we present an analytical approach for the form 
of P(o~) for low conductance regime. The distribution 
of the maximal random number a: max along the optimal 
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FIG. 3: Plots of P(a) for fixed (k)/a. The symbols are for 
the Kirchhoff method and the solid line is for the iterative 
method. For the same (k)/a, the iterative method for dif- 
ferent a shows the same P(cr) except that the lower cutoff is 
different, (a) ER network with (k)/a = 0.2. (b) ER network 
with (k)/a = 1.5. (c) SF network with (k)/a « 0.35, A = 2.5. 
The dashed line slopes are from the prediction of Eq. (|ll[) or 

m. 



path can be expressed in terms of the order parameter 
Pooip) m the percolation problem on the Cayley tree, 
where Pooip) is the probability that a randomly chosen 
node on the Cayley tree belongs to the IIPC [17]. For 
a random graph with degree distribution pf., the prob- 
ability to arrive at a node with k outgoing branche s by 
following a randomly chosen link is (k + l)pk/(k) [201 ] . 
The probability that starting at a randomly chosen link 
on a Cayley tree one can reach the £th generation is 



U(p) = u = i 



E 

fc=i 



p k k {\-pfi-x) 



fc-i 



(k) 



(4) 



where /o = 1. Slightly different from fi is the probability 
that starting at a randomly chosen node one can reach 
the nth generation, 



fn = 1 - ^Pk(l - Pfn-l) 



(5) 



fe=0 



In the asymptotic limit fi converges to for a given 
value of p, 



Pkk (I -pPoa) 



fc-1 



fe=l 



(k) 



(6) 



In this limit we have a pair of nodes on a random graph 
separated by a very long path of length n. The proba- 
bility that two nodes will be connected (conducting) at 
given p, can be approximated b y t he probability that 
both of them belong to the IIPC [II]: 



nip) 



PQ0(P) 

Poo(l) 



(7) 



where Poo(p) = hm^oo In = 1 - SfcloPfeC 1 ~ P p oo) k ■ 
Note that the negative derivative of II(p) with respect 
to p is the distribution of x max and thus gives P(cr) in 
the SD limit. In our case a — e~ ap , so replacing p by 
p = — hi a /a in Eq. ([7]) and differentiating with respect 
to a, we obtain the distribution of conductance in the SD 
limit when the source and sink are far apart (n — > oo), 



P(a) = --11(a) 



2Poo(p) dP^ip) 



aa[Poo(l)} 2 dp 



p=— In a /a 



(8) 

For ER networks the degree distribution is a Poisson 
distribution with p k = (k) k e -^ jk\ [li] and thusP oc (p) 
satisfies 

Poo(p) = l-e-^ p ^\ (9) 

which has a positive root Poo for p > p c = l/(k). And 
Poo(p) = Pooip), thus 



P(a) 



2P 00 (p) dPvoip) 



aalPUW dp 



\p— — In a I a i 



(10) 



where Poo(p) and Poo(l) are the solutions of Eq. ([9]). 

We test the analytical result Eq. (|10p by comparing 
the numerical solution of Eqs. §§§ and (jTUJ) with the sim- 
ulations on actual random graphs by solving Kirchhoff 
equations (Figs. [2] and [3]). The agreement between the 
simulations and the theoretical prediction is perfect in 
the SD limit, i.e. when (k)/a is small. 
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FIG. 4: Kirchhoff method results of the percolation part of 
ER networks with the same value of p c = l/{k) — 0.33. (a) 
Normalized P p (a) for fixed a/7V 1/3 = 1.5. (b) Scaled plot of 
(a)P(a) as function of a / (a) for three values of a/N 1 ^ 3 . For 
each value of a/N 1 ^ 3 , the thick line is for iV = 256 and the 
thin line is for N = 1024. 

Next we simplify P(cr) from Eq. (|10|) . Assuming that 
Poo(l) ~ 1 which is true for large (k) and approximating 
a slow varying function Poo(p) by P 00 (l) we obtain 



P(a) 



(11) 



for the range e~ a < a <C e~ apc with p c = l/(k). In 
Figs. [2] and [3] we also show the results predicted by 
Eq. ([TTj) . For an infinite network, for p < p c = 
Poo{p) = 0, and hence, the distribution of conductances 
must have a cutoff at a — e~ ap ". Indeed, in Fig.f^a) and 
Figs. [3£a) and(3Jb) we see that the upper cutoff of the 
iterative curves is close to e~ aPc . 
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As discussed above, the range of high conductivities 
corresponds to the case where both the source and the 
sink are on the IIPC. Previously we found this percola- 
tion part to scale as iV~ 2 / 3 . Using Fig. Ufa), we compute 
the integral for each P(a) from e~ aPc to oo, and find that 
indeed J^ aPc P(<r)dcr ~ N~ 2 / 3 , in good agreement with 
the theoretical approach. To show how the percolation 
part of P(cr) is related to the parameters N, a and p c , 
we analyze the conductance between pairs on the IIPC, 
i.e., each link on the optimal path from source to sink 
has x less than p c . We compute P p (<j) of these pairs on 
the IIPC. When we simulate this process, we have only 
N~ 2 /' 3 probability to find this part from the original nor- 
malized distribution P(cr). Thus, we normalize Pp(cr) by 
dividing by N~ 2 / 3 . Figures SJa) and|3Jb) show the nor- 
malized P p (cr) of pairs on the IIPC. In this range, we see 
that P p (cr) is dominated by high conductivities and we 
find (a) w e- a P-= and 

MW=/(^,^), (12) 

that is, for fixed apc/N 1 / 3 , (a)P p (a) scales with cr/(a) as 
seen in Fig.rjjb). The scaled distributions have the same 
shape for the same ap c /N x / 3 which specifies the strength 
of disorder similarly to the behavior of the optimal path 
lengths [H, QJ, [H, H3. The explanation of this fact for 
the distribution of conductances is analogous to the argu- 
ments presented in Refs. [13] and [IH for the distribution 



of the optimal path. Thus the position of the maximum 
of the scaled curves in Fig. [4jb) , and the whole shape of 
the distributions, depend on ap c /N 1 / 3 . 

We also find that the extreme high conductivities cor- 
respond to the case where source and sinks are separated 
by only one link. In this case, P(a) = ~ c -1 , 

(a<l). " 1 

In summary, we find that P(cr) exhibits two regimes. 
For a < e~ aPc , we show both analytically and numeri- 
cally that for ER networks P{a) follows a power law, 

P{cr)~<j- a [ a = 1- (k)/a\. (13) 

We also find that for SF networks, Eq. (fl~3|) seems to be 
a good approximation, consistent with numerical simula- 
tions. The distributions of optimal path length and the 
path length of the electrical currents in complex weighted 
networks [HI, [l!| have been found to depend on N for all 
length scales and all types of networks studied. In con- 
trast, here we find that the low conductance tail of P{<j) 
does not depend on N for both ER and SF networks. 
However, the high conductance regime (a > e~ apc ) of 
P(cr) does depend on N, in a similar way to the opti- 
mal path length and current path length distributions 

Qui. 
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